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The detection of 7-rays from the source HESS J1745— 290 in the Galactic Center (GC) region with 
the H.E.S.S. array of Cherenkov telescopes in 2004 is presented. After subtraction of the diffuse 
7-ray emission from the GC ridge, the source is compatible with a point-source with spatial extent 
less than 1.2'(stat.) (95% CL). The measured energy spectrum above 160 GeV is compatible with 
a power-law with photon index of 2.25 ± 0.04(stat.) ± 0.10(syst.) and no significant flux variation is 
detected. It is finally found that the bulk of the VHE emission must have non-dark-matter origin. 



PACS numbers: 98.70.Rz,98.35. Jk,95.35.+d 

INTRODUCTION 

Recently, the CANGAROO Q, VERITAS 0, 
H.E.S.S. H and MAGIC [i| collaborations have reported 
the detection of very high energy (VHE) 7-rays in the 
TeV energy range from the direction of the Galactic Cen- 



ter (GC). The nature of this source is still unknown. The 
main astrophysical explanations are based on particle ac- 
celeration in the region of the Sgr A East supernova rem- 
nant 0, in the vicinity of the supermassive black hole 
Sgr A* located at the center of our galaxy 0, 0] , or in a 
recently detected plerion Q. Another widely discussed 
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possibility concerns 7-ray emission from annihilation of 
dark matter (DM) particles [§]. 

Cosmological simulations of hierarchical structure for- 
mation 1(| ll[ predict that the DM particles form large 
scale structures in the Universe, and especially halos with 
a pronounced density cusp located at their center. Galax- 
ies are predicted to be embedded in such DM halos. Par- 
ticle physics and cosmology experiments constrain some 
characteristics of the new particles [IH : the new particles 
should be massive (> some GeV) and have weak inter- 
actions with ordinary matter of the order of the electro- 
weak cross sections. 

Extensions of the standard model of particle physics 
provide new particle candidates consistent with cosmo- 
logical DM and are of main interest to solve both is- 
sues. These models include supersymmetric theories (e.g. 
MSSMH [13 or AMSBHQ) or Kaluza-Klein (KK) 
scenarios with extra-dimensions [TH- 
AU DM particle candidates have some common proper- 
ties that can be used to detect them indirectly, since their 
annihilation may give rise to 7-rays, but also to neutrinos 
and cosmic-rays. Their annihilation rate is proportional 
to the square density of DM. It is thus enhanced in the 
dense DM regions at the center of DM halos. Cuspy ha- 
los may therefore provide detectable fluxes of VHE 7-rays 
(see 12j and references therein). The centers of galaxies 
are indeed good candidates for indirect DM detection, 
the closest candidate being the center of the Milky Way. 
The 7-ray energy spectrum generated by DM annihila- 
tion is characterized by a continuum ranging up to the 
mass of the DM particle, and possibly faint 7-ray lines 
provided by two-body final states 0, [l]| ■ 

For annihilation of DM particles of mass todm accumu- 
lated in a spherical halo of mass density profile p(r) and 
particle density profile p(r)/moM, the 7-ray flux F{E) is 
proportional to the line-of-sight-integrated squared par- 
ticle density, multiplied by the velocity-weighted anni- 
hilation cross section (av) and the number of photons 
dN^/dE generated per annihilation event can 
be factored into a term J depending on the halo parame- 
ters and a term depending on the particle physics model: 
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with F = 2.8 • 10 _12 cni _2 s _1 . J(Afi) denotes the ave- 
rage of J over the solid angle Af2 corresponding to the 
angular resolution of the instrument, normalized to the 
local DM density 0.3 GeV cm" 3 : 
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The shape of the measured 7-ray spectrum depends only 
on the particle properties, embedded in the term F(E), 
and especially on the 7-ray multiplicity dN^/dE. The 
measured angular distribution of the 7-rays depends only 



on J(Af2). The overall 7-ray flux depends on both terms. 
Close to the GC, halo density profiles are predicted to fol- 
low a power-law ph(t) ~ r~ a with a between 1 [loj and 
1.5 Recent N-body simulations 17| suggest that a 
could monotonically decrease to zero towards the GC. 
Values of a lower than ~ 1.2 lead to an angular distri- 
bution broader than the H.E.S.S. angular resolution and 
can thus be constrained. 

In this letter, we present results on VHE 7-rays from 
the GC based on a dataset collected in 2004 with the com- 
plete H.E.S.S. array of imaging atmospheric Cherenkov 
telescopes (IACTs). 



THE H.E.S.S. TELESCOPES AND THE 
GALACTIC CENTER DATA SET 

H.E.S.S. is a system of four IACTs (see [5 and ref- 
erences therein) located in Namibia, close to the Tropic 
of Capricorn. The telescopes stand at the corners of a 
square of 120 m side. The Cherenkov light emitted by 
7-induced air showers is imaged onto cameras of 960 pho- 
tomultipliers, covering a field of view of 5° in diameter. 
The large mirror area of 107 m 2 per telescope results in 
an energy threshold of 100 GeV at Zenith The 
stereoscopic imaging of air showers allows the precise re- 
construction of the direction and energy of the 7-rays. 

The previously published H.E.S.S. results on the GC Q 
were based on 17 h of data recorded with the first two 
telescopes in 2003. Here, we report on results obtained 
with the full four-telescope array, using 48.7 h (live time) 
of data collected between March 30th and September 4th, 
2004. The full array provides higher detection rates than 
the 2003 data, as well as improved background rejection 
and angular resolution. The bulk of the data (33.5 h) 
were obtained in "wobble mode" , where the source region 
is displaced by typically ±0.7° from the optical axis of 
the system. An additional 15.2 h dataset was obtained 
from the Galactic plane survey [l||, within 2° of Sgr A*. 

Two different techniques for calibration and image 



analysis were applied 1 



20] and give identical results. 
Both methods provide a typical energy resolution of 15% 
and an angular resolution of 0.1° above the analysis en- 
ergy threshold. The results described in this paper are 
derived using the second technique. 



TEV 7-RAYS FROM THE DIRECTION OF THE 
GALACTIC CENTER 



The data show an excess of 1863 7-events from 
HESS J1745-290 within 0.1° from the GC (see 
Fig. Q] )• This excess is detected on top of a 
hadronic background of 1698 events, with a significance 
of 37.9 standard deviation above background, calculated 
according to [21 1. 
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Diffuse 7-ray emission extended along the galactic 
plane has been discovered in these data and was reported 
elsewhere [22]. It was shown that this emission likely 
originates in cosmic-ray interactions with giant molecular 
clouds and is thus proportional to the density of cosmic- 
rays and of target material. To study the shape and po- 
sition of HESS J1745— 290, the diffuse emission has been 
modeled assuming a perfect correlation with the molecu- 
lar cloud density from CS data ji^. Cosmic-ray density 
was assumed to have a Gaussian dependence on distance 
to the GC with scale a — 0.8°. The resulting emission 
model has been smeared with the H.E.S.S. PSF (point 
spread function, approximately Gaussian with a 68% 
containment radius of 0.1°). The H.E.S.S. central source 
has been fitted as a superposition of the diffuse compo- 
nent and either a point-like source, a Gaussian source or a 
DM halo shape. Likelihood fits of these different models 
to the 7-ray count-map within a radius of 0.5° of Sgr A* 
were made with the flux normalisation of the diffuse emis- 
sion model as a free parameter. Assuming a point source 
for HESSJ1745-290, folded with the H.E.S.S. PSF, the 
best fit location of the source is (£ = 359°56'33.3" ± 
9.7", b = -0°2'40.6" ± 10") in Galactic coordinates 
or (a = 17 h 45 m 39.44 s ± 0.6 s , 6 = -29d00'30.3" ± 9.7") 
in equatorial coordinates (J2000.0), within 7" ± 14" tat ± 
28" yst from the putative supermassive black hole Sgr A* . 
Improvements in the pointing accuracy may allow the 
systematic errors to be reduced in the future. No remain- 
ing contribution is found in the 7-ray map after subtrac- 
tion of the fitted emission, indicating that this model is 
consistent with the data. The distribution of the angle 8 
between the 7-ray direction and the position of Sgr A* 
after subtraction of the fitted diffuse emission is shown 
in Fig. [1] and is consistent with the H.E.S.S. PSF. The 
diffuse emission is found to contribute to 16% of the to- 
tal signal of HESS J1745-290 within 0.1°. Assuming a 
azimuthally symmetric Gaussian brightness distribution 
centered on the best fit position given above, folded by 
the H.E.S.S. PSF, an upper limit on the source size of 1.2' 
(95% CL) was derived (including statistical errors only). 

The compatibility of the spatial extension of 
HESS J1745-290 with a DM halo centered on Sgr A* 
and with density following p(r) ex r~ Q was also tested. 
Different values of the logarithmic slope a were assumed. 
The diffuse component and the DM halo were both folded 
with the H.E.S.S. PSF. Leaving both normalisations free, 
the fit likelihood is compared to the point-like source hy- 
pothesis discussed above in order to derive a lower limit 
on the slope a of 1.2 (95% CL). 

The spectral energy distribution (SED) of 7-rays F(E) 
of the GC source is determined using an 0.1° integra- 
tion radius and assuming a point source. As the flux 
contamination of the diffuse emission (16%) is of the 
same order as flux systematic errors, it is not subtracted 
in this analysis. Moreover, as the shape of the diffuse 
emission spectrum is compatible with that of the central 
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FIG. 1: (Color online) Background-subtracted distribution 
of the angle 9 between the 7-ray direction and the position 
of Sgr A*. Circles: all detected 7-rays events. Open tri- 
angles: central object after subtraction of the 7-ray diffuse 
emission model (see text). Line: calculated PSF normalized 
to the number of 7-rays within 0.1° after subtraction is also 
shown. The distribution of events after subtraction matches 
the calculated PSF while the initial distribution shows a sig- 
nificant tail. The variation of the PSF related to the source 
energy spectrum, zenith angle and offset position in the field 
of view are taken into account. Insert: same distribution for 
the point-like source PKS 2155-304 0. The calculated PSF 
(line) also matches the data. 
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FIG. 2: (Color online) Spectral energy density E 2 x (IN /(IE 
of 7-rays from the GC source, for the 2004 data (full points) 
and 2003 data || (open points). Upper limits are 95% CL. 
The shaded area shows the power-law fit dN/dE ~ E~ T . 
The dashed line illustrates typical spectra of phenomenologi- 
cal MSSM DM annihilation for best fit neutralino masses of 
14 TeV. The dotted line shows the distribution predicted for 
KK DM with a mass of 5 TeV. The solid line gives the spec- 
trum of a 10 TeV DM particle annihilating into t + t~ (30%) 
and bb (70%). 



source 22], the measured spectral shape is not altered. 
The SED is shown in Fig. [2] (together with the spec- 
trum derived from the H.E.S.S. 2003 data). Although 
a 7-ray excess is seen at energies as low as 100 GeV, 
the spectrum shown is calculated only above 160 GeV 
to eliminate systematic errors arising from an energy 
reconstruction bias close to threshold. Over the en- 
ergy range 160 GeV — 30 TeV the energy spectrum can 
be characterized by a power-law, F(E) ~ E~ T with 
T = 2.25 ± 0.04(stat.) ± 0.10(syst.) (with a fit proba- 
bility of 39%). The 2003 and 2004 spectra are consistent 
in shape and normalization, with an integral flux above 
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1 TeV of (1.87±0.10(stat.)±0.30(syst.)) x 1CT cm-'s" 
There is no evidence for a cut-off in the spectrum and 
lower limits at 95% CL of 9 TeV and 7 TeV are derived 
assuming an exponential cut-off and a sharp cut-off (38[, 
respectively. The experimental spectrum has also been 
fitted as a sum of a free power-law and a monoenergetic 
7-ray line [39| whose energy and normalisation have been 
scanned. No indications of line emission are found. 

There is no significant variation in flux between 2003 
and 2004; data are consistent with a constant flux [25j ]. 
Searches for variability or flares on time scales down to 
10 min did not show statistically significant deviations 
from the mean flux. We note that approximately 20 min 
of data are required for a 3 standard deviation detec- 
tion of the source above background. A flare lasting for 
10 min (30 min, 3 h, respectively) and with a 7-fold (4- 
fold, 2-fold, respectively) increase over the quiescent flux 
would be detected at the 99% CL. Data were also ana- 
lyzed for periodic or quasi-periodic variations on scales 
between 1 mHz and 16 /HIz, using the Lomb-Scargle 
method [2^|. Again, no statistically significant period- 
icity was found. However, if the VHE emission is asso- 
ciated with Sgr A* and given its typical rate of X-ray 
flares of 1.2 per 24 h {27], the 48.7 h of H.E.S.S. data may 
simply not contain a flare event. 



DARK-MATTER INTERPRETATION 

The location of the TeV 7-ray signal and its tempo- 
ral stability are consistent with a DM annihilation signal 
from a halo centered Sgr A*. 

In a first step, it is assumed that all 7-rays from 
HESS J1745-290 are due to DM annihilations. The hy- 
pothetical DM halo centered on Sgr A* was found in 
the previous section to be very cuspy, with a logarithmic 
slope a higher than 1.2. This value is consistent with 
the DM halo shapes predicted by some structure forma- 
tion simulations. The energy spectrum provides another 
crucial test concerning a possible DM origin for the de- 
tected VHE emission. The extension of the spectrum 
beyond 10 TeV requires masses of DM particles which 
are uncomfortably large MSSM. The annihilation spec- 
tra of phenomenological MSSM neutralinos depend on 
the gaugino/higgsino mixing, but all exhibit a curved 
spectrum, which in a E 2 dN/dE representation rises for 
E <C ttidm, plateaus at E/tubm ~ 0.01 — 0.1, and falls off 
approaching ttidm- AMSB models lead to similar spec- 
tra. Such a spectral shape is inconsistent with the mea- 
sured power-law as seen in Fig. [51 H.E.S.S. data from 
2003, with restricted energy range and lower statistics, 
were still marginally consistent with DM spectra [28| , but 
it appears impossible to generate a power-law extending 
over two decades from the quark and gluon fragmentation 
spectra of neutralino decays, also considering radiative 
effects 1291. As an alternative scenario, mixed t + t~, bb 



final states have been proposed [30(, with DM masses in 
the 6 — 30 TeV range, generating a flatter spectrum. Non- 
minimal SUSY models can be constructed which allow 
such decay branching ratios, combined with neutralino 



masses of tens of TeV. KK DM discussed in [3l| also 
give harder spectra. PYTHIA 6.225[32j was used to com- 
pute the contributions from all annihilation channels [40l ] . 
However, all the tested model spectra still deviate signif- 
icantly from the observed power-law spectrum as shown 
in Fig. H 

On the other hand, if the bulk of the VHE emission 
has non-DM origin, there is still the possibility of a DM 
signal hidden under an astrophysical spectrum. To search 
for such a contribution, we fitted the experimental 
spectrum as the sum of a power-law with free normal- 
ization and index, and a MSSM (or KK) spectrum. 
Leaving the normalisation of the DM signal free, the 
range of ttidm is scanned. For the MSSM, annihilation 
spectra dN^/dE = No/m x (E/m x ) exp ( — cE/m x ^ 
are used with three different sets of parameters, 
one approximating the average annihilation spec- 
trum ((N ,T,c) = (0.081,2.31,4.88)) the other two 
((iVo,r,c) = (0.2,1.7,10) and (0.4,1.7,3.5)) roughly 
encompassing the range of model spectra generated 



using Dark Susy [16| [4lJ. No significant DM component 



is detected with this procedure, the DM component flux 
upper limit being of the order of 10% of the source flux. 
Assuming a NFW profile, 99% CL upper limits on the 
velocity-weighted annihilation cross section (av) are of 
the order of 10~ 24 — 10 _23 cm 3 s^ 1 , above the predicted 
values of the order of 3 x 10~ 26 cm 3 s _1 . These limits can 
vary by plus or minus three orders of magnitude if one 
assumes other DM halo shapes. In the case of adiabatic 
compression of DM due to the infall of baryons to the 
GC, the flux could be boosted up to a factor 1 000 [33j |. 
The H.E.S.S. data might then start to exclude some (av) 
values. 

In conclusion, the power-law energy spectrum of the 
source HESS J1745-290 measured using the H.E.S.S. tele- 
scopes show that the observed VHE 7-ray emission is not 
compatible with the most conventional DM particle an- 
nihilation scenarios. It is thus likely that the bulk of the 
emission is provided by astrophysical non-DM processes. 
However, due to high density of candidate objects for 
non-thermal emission within the source region the na- 
ture of the source is not clear. 



* Also at Institut fur Physik, Humboldt-Universitat zu 
Berlin, Newtonstr. 15, D 12489 Berlin, Germany 

t Also at European Associated Laboratory for Gamma- 
Ray Astronomy, jointly supported by CNRS and MPG 

* Also at Landessternwarte, Konigstuhl, D 69117 Heidel- 
berg, Germany 



5 



§ Also at Institut fur Astronomie und Astrophysik, Uni- 

versitat Tubingen, Sand 1, 72076 Tubingen, Germany 
" Also at Max-Planck-Institut fiir Kernphysik, P.O. Box 

103980, D 69029 Heidelberg, Germany 
** Now at Purdue University, Department of Physics, 525 

Northwestern Avenue, West Lafayette, IN 47907-2036, 

USA 

^ Electronic address: ripkenj@mail.desy.de 

** Electronic address: rollandl@in2p3.fr; Also at DAP- 

NIA/DSM/CEA, CE Saclay, F-91191 Gif-sur-Yvette, 

Cedex, France 

§§ URL: |http : //www . mpi-hd ■ mpg . deTh f m/HESS/HESS . html] 
[1] K. Tsuchiya et al. (CANGAROO Collaboration), ApJ 

606, L115 (2004). 
[2] K. Kosack et al. (VERITAS Collaboration), ApJ 608, 

L97 (2004). 

[3] F. Aharonian et al.(HESS Collaboration), A&A 425, L13 

(2004) . 

[4] J. Albert et al. (MAGIC Collaboration), ApJ 638, L101 
(2006). 

[5] R.M. Crocker et al., ApJ 622, 892 (2005). 
[6] F. Aharonian & A. Neronov, ApJ 619, 306 (2005). 
[7] A. Atoyan & CD. Dermer, ApJ 617, L123 (2004). 
[8] Q.D. Wang, MNRAS 367, 937 (2006). 
[9] L. Bergstrom, Rep. Prog. Phys. 63, 793 (2000). 
[10] J.F. Navarro, C.S. Frenk & S.D.M. White, ApJ 490, 493 
(1997). 

[11] B. Moore et al., MNRAS 310, 1147 (1999). 

[12] G. Bertone, D. Hooper & J. Silk, Phys. Rep. 405, 279 

(2005) . 

[13] J. Ellis et al., Eur. Phys. J. C24, 311 (2002). 

[14] S. Profumo & P. Ullio, JCAP 07, 006 (2004). 

[15] G. Servant & T.M. Tait, Nucl. Phys. B650, 391 (2003). 

[16] P. Gondolo et al., JCAP 07, 008 (2004). 

[17] J.F. Navarro et al., MNRAS 355, 794 (2004). 



[18] F. Aharonian et al., Astron. Astrophys. 457, 899 (2006). 
[19] F. Aharonian et al., Science 307, 1938 (2005). 
[20] L. Rolland & M. de Naurois, AIP Conf. Proc. 745, 715 
(2004). 

[21] T. Li & Y. Ma, ApJ 272, 317 (1983). 
[22] F. Aharonian et al., Nature 439, 695 (2006). 
[23] M. Tsuboi, et al., ApJS 120, 1 (1999). 
[24] F. Aharonian et al., Astron. Astrophys. 442, 895 (2005). 
[25] L. Rolland et al., Proc. 29th ICRC, Pune (2005). 
[26] J.D. Scargle, ApJ 343, 874 (1989). 
[27] M.P. Muno et al., ApJ 589, 225 (2003). 
[28] D. Horns, Phys. Lett. B607, 225 (2005). 
[29] L. Bergstrom et al., Phys. Rev. Lett. 95, 241301 (2005). 
[30] S. Profumo, Phys. Rev. D72, 10352 (2005). 
[31] L. Bergstrom et al., Phys. Rev. Lett. 94, 131301 (2005). 
[32] T. Sjostrand et al., Computer Phys. Commun. 135, 238 
(2001). 

[33] F. Prada et al., Phys. Rev. Lett. 93, 241301 (2004). 

[34] MSSM: Minimal Supersymmetric Standard Model. 

[35] AMSB: Anomaly Mediated Supersymmetry Breaking. 

[36] The energy threshold is defined as the peak of the differ- 
ential reconstructed energy distribution for a 7-ray source 
with a power-law energy spectrum with photon index 2.6. 

[37] From a dN/d9 excess distribution, the y- values have been 
divided by 29 (9 being the bin center). 

[38] dN/dE = above the cut-off energy. 

[39] The monoenergetic 7-ray line shape is estimated as the 
H.E.S.S. energy resolution. Its variation related to the 
line-energy, zenith angle, and position of the source in 
the field of view, are taken into account. 

[40] Comparison of different versions of PYTHIA has shown 
that there are systematic uncertainties of the order of 
10% in the fluxes and the spectral shapes. 

[41] DarkSUSY version 4.1. 



